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1  Introduction 

The  aim  of  the  report  is  to  show  novel  methods  in  functional  neuroimaging  [1] 
of  sleep  inertia  stages  with  emphasis  on  novel  time-frequency  preprocessing  meth- 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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ods  [2,  3,  4,  5,  6]  and  models  of  effective  connectivity  evaluated  by  dynamic  causal 
modeling  and  econometric  models  (such  as  Granger  causality  and  related  meth¬ 
ods)  [7,  8,  1] 

2  Methods 

The  EEG  signals  from  two  conditions  4^1  and  for  all  experimental  factors 
EEGQS  till  EEG17  for  the  following  channels:  G4  —  Al;  G3  —  A2;  Ez  —  Al] 
Pz  —  Al;  Oz  —  A2  where  hrst  preprocessed  with  EMD  adaptive  hltering  procedures 
(removing  mostly  muscle  artifacts  in  lower  and  high  frequency  ranges).  The  so 
obtained  clean  EEG  signals  were  further  with  interdependence  measures. 

2.1  Empirical  Mode  Decomposition  (EMD)  EEG  Prepro¬ 
cessing 

...to  be  added  soon... 

2.2  Signals  Interdependence  Evaluation 

Suppose  that  we  are  given  n  signals  Xi{k),  X2{k), . . . ,  Xn{k),  each  stemming  from 
a  different  channel.  We  consider  the  multivariate  autoregressive  (MVAR)  model; 

.Y(l-)  =  ^A(j)A'(i-i)  +  £W,  (1) 

1=1 

where  X{k)  =  {Xi{k),  X2{k), . . .  ,Xn{k))^,  p  is  the  model  order,  the  model  coeffi¬ 
cients  A(j)  are  n  x  n  matrices,  and  E{k)  is  a  zero-mean  Gaussian  random  vector 
of  size  n.  In  words:  Each  signal  Xi{k)  is  assumed  to  linearly  depend  on  its  own 
p  past  values  and  the  p  past  values  of  the  other  signals  Xj{k).  The  deviation  be¬ 
tween  X{k)  and  this  linear  dependence  is  modeled  by  the  noise  component  E{k). 
Model  (1)  can  also  be  cast  in  the  form: 

BW  =  EA(j)A(fc-;),  (2) 

1=0 

where  A(0)  =  I  (identity  matrix)  and  A(j)  =  — A(j)  for  j  >  0.  One  can  transform 
(2)  into  the  frequency  domain  (by  applying  the  ;2-transform  and  by  substituting 
A  where  1/At  is  the  sampling  rate): 

A(/)=A-i(/)E(/)  Ah(/)E(/).  (3) 
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The  power  spectrum  matrix  of  the  signal  X{k)  is  determined  as 

S(/)  ^  X{f)X{fr  =  H(/)VH(/),  (4) 

where  V  stands  for  the  covariance  matrix  of  E{k). 


2.2.1  Coherence  (COH) 

The  coherence  function  quantifies  linear  correlations  in  frequency  domain.  One 
distinguishes  the  magnitude  square  coherence  function  and  the  phase  coherence 
fnnction.  The  former  is  dened  as: 

,  \X{f)Y*{f)f 

\x{f)\\Yur 

where  X{f)  and  Y (/)  are  the  Fourier  transforms  of  x  and  y  respectively;  Y*  is  the 
complex  conjugate  oiY  G  C,  and  |1^|  is  the  magnitnde  of  Y  .  The  phase  coherence 
fnnction  is  dened  as 

<#.(/)  ^arg[X(/)y(/)].  (6) 

In  practice,  one  often  snbdivides  the  signals  x  and  y  in  M  segments  (of  eqnal 
length),  and  determines  c(/)  by  averaging  over  those  segments: 

,  \{X{f)Y*{f))\^ 

mm\{Y{f))y 

where  (■)  denotes  averaging  over  the  M  segments.  Along  the  same  lines,  the  phase 
coherence  0(/)  is  often  compnted  as: 

m  A  arg|{A'(/)y(/)>|.  (8) 


(5) 


Note  that  both  c(/)  and  0(/)  depend  on  the  frequency  /. 

As  for  properties  of  coherence,  its  estimated  value  ranges  between  0  and  1.  For 
a  given  freqnency.  A  valne  of  0  indicates  that  the  activities  of  the  signals  in  this 
particular  frequency  bin  are  linearly  independent,  whereas  a  maximnm  value  of  1 
gives  the  top  linear  correlation  for  snch  particular  frequency  bin. 


2.2.2  Directed  Transfer  Function  (DTF) 

\H^m" 


^2, 

hj 


(/)  = 


Em 
7  =  1 


^q(/)r 


e  [0,1], 


(9) 


where  the  (frequency-dependent)  normalization  is  chosen  so  that  'yfj{f)  quanties 
the  fraction  of  inow  to  channel  i  stemming  from  channel  j. 
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(10) 


2.2.3  Partial  Directed  Coherence  (PDC) 


PM  = 


Mf) 


Mf) 


2 


e  C, 


2.2.4  Direct  Directed  Transfer  Function  (dDTF) 

X?,(/)  =  f)?(/)Cj(/)€lO,l],  (11) 

which  is  non-zero  if  the  connection  between  channel  i  and  j  is  causal  (non-zero 
Pijif))  direct  (non-zero  Cfj{f)). 


3  Preliminary  pairwise  channel  synchrony  eval¬ 
uation  results 

3.1  Magnitude  coherence 
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Figure  1:  Magnitude  coherence  as  in  (5):  Condition  &  factor  EEG08.  X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  2:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG08.  X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
hltered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  3:  Magnitude  coherence  as  in  (5):  Condition  &  factor  EEG09.  X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  4:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG09. X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  5:  Magnitude  coherence  as  in  (5):  Condition  &  factor  EEG10.X-aa:e5 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  6:  Magnitude  coherence  as  in  (5):  Condition  ^^3  &  factor  EEG10.X-tta:e5 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 


7 


C4-A1 


1 0  20  30  40 
C3-A2 


1 

0.5 

0 


1 0  20  30  40 
Pz-AI 


1 0  20  30  40 
Fz-A1 


1 0  20  30  40 
EEG11:1:  rCOH 


Figure  7:  Magnitude  coherence  as  in  (5):  Condition  7^1  &  factor  EEGll.X-aa:e5 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  8:  Magnitude  coherence  as  in  (5):  Condition  ^^3  &  factor  EEGll.X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  9:  Magnitude  coherence  as  in  (5):  as  in  (5):  Condition  &  factor 
EEG12. X-axes  in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around 
IHz  are  usually  filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  10:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG12. X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
hltered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  11;  Magnitude  coherence  as  in  (5):  Condition  ^^1  &  factor  FjEGlS.X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 


C3-A2 


10  20  30  40 
Fz-A1 


C4-A1  C4-A1 


1  I - 1  1 


1 0  20  30  40 


Fz-AI 


C4-A1  C4-A1 


10  20  30  40  10  20  30  40 

C3-A2  C3-A2 


1  I - 1  1 


10  20  30  40  10  20  30  40 

Fz-AI  Fz-AI 


_ :  qI _ I  qI _ I  --cmm 

10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 

Pz-AI  Pz-AI  Pz-AI  Pz-AI 


< 


O 


< 


o 


10  20  30  40 
Oz-A2 


10  20  30  40 


0^ - O' - 1 

10  20  30  40  10  20  30  40 

Oz-A2  Oz-A2 


1| - 1  1 


10  20  30  40  10  20  30  40 


1 

CJ 


1 0  20  30  40 


Oz-A2  EEG13:3;rCOH 


10  20  30  40 


Figure  12:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG13. X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  13:  Magnitude  coherence  as  in  (5):  Condition  &  factor  EEGIA. X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  14:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG14. X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  15:  Magnitude  coherence  as  in  (5):  Condition  &  factor  EEG15. X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  16:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG15. X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
hltered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  17:  Magnitude  coherence  as  in  (5):  Condition  &  factor  FjEGlG.X-axes 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  18:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG16. X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  19:  Magnitude  coherence  as  in  (5):  Condition  &  factor  EEG17.X-aa:e5 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  20:  Magnitude  coherence  as  in  (5):  Condition  ^3  &  factor  EEG17.X-aa;es 
in  all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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3.2  Phase  coherence 
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Figure  21:  Phase  coherence  as  in  (6):  Condition  &  factor  EEG08.  X-axes  in 
all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  22:  Phase  coherence  as  in  (6):  Condition  ^^3  &  factor  EEG08.  X-axes  in 
all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  23:  Phase  coherence  as  in  (6):  Condition  &  factor  EEG09.  X-axes  in 
all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  24:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG09.  X-axes  in 
all  plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually 
filtered  by  EMD  thus  no  activity  there  analyzed). 
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Figure  25:  Phase  coherence  as  in  (6):  Condition  7^1  &  factor  FjFjGIO. X-axes  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  26:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEGIO. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  27:  Phase  coherence  as  in  (6):  Condition  7^1  &  factor  EEGll.X-aa:e5  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  28:  Phase  coherence  as  in  (6):  Condition  ^^3  &  factor  EEGll. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  29:  Phase  coherence  as  in  (6):  Condition  &  factor  FjFjG12.X-axes  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  30:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG12. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  31:  Phase  coherence  as  in  (6):  Condition  7^1  &  factor  FjEGlS.X-axes  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  32:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG13. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  33:  Phase  coherence  as  in  (6):  Condition  &  factor  EEG14.X-aa:e5  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  34:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG14. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  35:  Phase  coherence  as  in  (6):  Condition  &  factor  FjFjG15.X-axes  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  36:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG15. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  37:  Phase  coherence  as  in  (6):  Condition  &  factor  EEGIG. X-axes  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  38:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG16. X-axes  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  39:  Phase  coherence  as  in  (6):  Condition  &  factor  EEG17.X-aa;es  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered 
by  EMD  thus  no  activity  there  analyzed). 
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Figure  40:  Phase  coherence  as  in  (6):  Condition  ^3  &  factor  EEG17.X-aa;es  in  all 
plots  are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered 
by  EMD  thus  no  activity  there  analyzed). 
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3.3  Directed  transfer  function  (DTF) 


26 


10  20  30  40 
C3-A2 


10  20  30  40 
C3-A2 


10  20  30  40 
C3-A2 


10  20  30  40 
Fz-AI 


10  20  30  40 
Pz-AI 


10  20  30  40 
Pz-AI 


10  20  30  40 
Fz-AI 


10  20  30  40 
Fz-AI 


10  20  30  40 
Pz-AI 


10  20  30  40 
Oz-A2 


10  20  30  40  10  20  30  40 

Oz-A2  Oz-A2 

- 1  1  I - 1 


10  20  30  40 
EEG08:1:  DTF 


Figure  41:  DTF  as  in  (9):  Condition  and  factor  EEG08.  X-axes  in  all  plots 
are  frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by 
EMD  thus  no  activity  there  analyzed). 
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Figure  42:  DTF  as  in  (9):  Condition  ^3  and  EEG08.  X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  43:  DTF  as  in  (9):  Condition  and  factor  EEG09.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  44;  DTF  as  in  (9):  Condition  ^^3  and  EEG09. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  45:  DTF  as  in  (9):  Condition  and  factor  EEG10.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  46;  DTF  as  in  (9):  Condition  ^3  and  EEGIO. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  47:  DTF  as  in  (9):  Condition  and  factor  FjEGll.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  48;  DTF  as  in  (9):  Condition  ^3  and  EEGll. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  49:  DTF  as  in  (9):  Condition  and  factor  FjFjG12.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  50:  DTF  as  in  (9):  Condition  ^3  and  EEG12. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  51:  DTF  as  in  (9):  Condition  and  factor  EEGIS.X- axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  52;  DTF  as  in  (9):  Condition  ^3  and  FjEGlS.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  53:  DTF  as  in  (9):  Condition  and  factor  EEGl^.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  54;  DTF  as  in  (9):  Condition  ^3  and  EEG14. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  55:  DTF  as  in  (9):  Condition  and  factor  YEjGlb.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  56;  DTF  as  in  (9):  Condition  ^3  and  EEG15. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  57:  DTF  as  in  (9):  Condition  and  factor  YEjGlQ.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  58:  DTF  as  in  (9):  Condition  ^^3  and  EEG16.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 


35 


C4-A1  C4-A1  C4-A1 


10  20  30  40  10  20  30  40  10  20  30  40 

C3-A2  C3-A2  C3-A2 


10  20  30  40  10  20  30  40  10  20  30  40 

Oz-A2  Oz-A2  Oz-A2  Oz-A2 


10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


1 

0.5 
0 

1 

0.5 
0 

1 

0.5 
0 

1 

0.5 
0 

10  20  30  40 
EEG17:1:  DTF 


C4-A1 


Figure  59:  DTF  as  in  (9):  Condition  and  factor  YEjGVJ .X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  60;  DTF  as  in  (9):  Condition  ^^3  and  Yj'EGVJ .X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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3.4  Direct  directed  transfer  function  (dDTF) 
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Figure  61:  dDTF  as  in  (11):  as  in  (9):  Condition  ^^1  and  EEG08.  X-axes  in  all 
plots  are  frequencies  0  — 20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered 
by  EMD  thus  no  activity  there  analyzed). 


10203040 

C3-A2 


10203040 

C3-A2 


10203040 

C3-A2 


10203040 

Fz-AI 


10203040 

Fz-A1 


10203040 

Fz-AI 


10203040 

Pz-A1 


10203040 

Pz-A1 


10203040 

Pz-AI 


10203040 

Oz-A2 


10203040 

Oz-A2 


10203040 

Oz-A2 


10203040 
EEG08;3;  dDTF 


Figure  62:  dDTF  as  in  (11):  Condition  ^(13  and  EEG08.  X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  63:  dDTF  as  in  (11):  Condition  ^^1  and  EEG09.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  64:  dDTF  as  in  (11):  Condition  ^(13  and  EEG09. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  65:  dDTF  as  in  (11):  Condition  ^^1  and  FjFjGIO. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  66:  dDTF  as  in  (11):  Condition  ^3  and  FjFjGIO. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  67:  dDTF  as  in  (11):  Condition  ^^1  and  FjEGll.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  68:  dDTF  as  in  (11):  Condition  ^3  and  FjEGll.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  69:  dDTF  as  in  (11):  Condition  and  EEG12. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  70:  dDTF  as  in  (11):  Condition  7)^3  and  EEG12. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  71:  dDTF  as  in  (11):  Condition  ^^1  and  FjEGlS.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  72:  dDTF  as  in  (11):  Condition  ^3  and  EEG13. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  73:  dDTF  as  in  (11):  Condition  ^^1  and  EEGl^.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  74:  dDTF  as  in  (11):  Condition  ^3  and  EEGl^.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  75:  dDTF  as  in  (11):  Condition  and  FjEGlb.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  76:  dDTF  as  in  (11):  Condition  ^3  and  EEC  1 5. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  77:  dDTF  as  in  (11):  Condition  and  EEG16. X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  78:  dDTF  as  in  (11):  Condition  ^3  and  EEC  1 6. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  79:  dDTF  as  in  (11):  Condition  ^^1  and  EEG17.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  80:  dDTF  as  in  (11):  Condition  ^^3  and  YEjGII .X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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3.5  Partial  directed  coherence  (PDC) 
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Figure  81:  PDC  as  in  (10):  Condition  and  EEG08.  X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  82:  PDC  as  in  (10):  Condition  ^3  and  EEG08.  X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  83:  PDC  as  in  (10):  Condition  and  EEG09.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  84:  PDC  as  in  (10):  Condition  ^(13  and  EEG09. X-aa;es  in  all  plots  are 
freqnencies  0  —  20Hz  (DC  and  freqnencies  aronnd  IHz  are  usnally  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  85:  PDC  as  in  (10):  Condition  and  FjEGlO.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  86:  PDC  as  in  (10):  Condition  ^3  and  EEGIO. X-aa;es  in  all  plots  are 
freqnencies  0  —  20Hz  (DC  and  freqnencies  aronnd  IHz  are  usnally  filtered  by  EMD 
thus  no  activity  there  analyzed). 


51 


C4-A1 


C4-A1 


C4-A1 


C4-A1 


C3-A2 


1 


10  20  30  40 


1 


10  20  30  40 

Fz-AI  Fz-A1 


CJ 

pifiiiHii 

N 

Q. 

0 

N 

O 

0 

ilHlIii 

10  20  30  40 
Fz-AI 


10  20  30  40 
Fz-AI 


1| - 1  1 


10  20  30  40  10  20  30  40 

Pz-A1  Pz-AI  Pz-AI 


BWilllfflifflfflnf 

C\J 

< 

CO 

o 

0.5 

1  0-5 

LL 

. !"ii4mimiiiiiiiiiiiirttt 

0 

0 

[tftfffftppiniiiiiiiiiiiiiiiH# 

10  20  30  40  10  20  30  40  10  20  30  40 

Oz-A2  Oz-A2  Oz-A2 


C\J 

< 

< 

< 

cJ,  0.5 

^  0-5 

o 

o 

LL 

0 

0 

0 

10  20  30  40  10  20  30  40  10  20  30  40 


1 


10  20  30  40 


Oz-A2 


10  20  30  40 


1 


0.5 


10  20  30  40 
Pz-AI 


1 


0.5 


10  20  30  40 
EEG11:1;  PDC 


Figure  87:  PDC  as  in  (10):  Condition  and  FjEGll.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  88:  PDC  as  in  (10):  Condition  ^3  and  FjEGll.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  89:  PDC  as  in  (10):  Condition  and  EEG12. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  90:  PDC  as  in  (10):  Condition  ^3  and  EEG12. X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  hltered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  91:  PDC  as  in  (10):  Condition  ^^1  and  EEGIS.X- axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  92:  PDC  as  in  (10):  Condition  ^3  and  FjEGlS.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 


54 


C4-A1  C4-A1  C4-A1  C4-A1 


10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


C3-A2  C3-A2  C3-A2  C3-A2 


10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


Oz-A2  Oz-A2  Oz-A2  Oz-A2  EEG14:1;PDC 


10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


Figure  93:  PDC  as  in  (10):  Condition  and  EEG14.X-aa:e5  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 


C4-A1  C4-A1  C4-A1 


10  20  30  40  10  20  30  40  10  20  30  40 


C3-A2  C3-A2  C3-A2 


10  20  30  40  10  20  30  40  10  20  30  40 

Pz-A1  Pz-AI  Pz-AI 


10  20  30  40  10  20  30  40  10  20  30  40 

Oz-A2  Oz-A2  Oz-A2  Oz-A2 


10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


C4-A1 


10  20  30  40 
C3-A2 


10  20  30  40 
Fz-AI 


10  20  30  40 
Pz-AI 


10  20  30  40 
EEG14;3;  PDC 


Figure  94:  PDC  as  in  (10):  Condition  ^^3  and  Yj'EGlA.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  95:  PDC  as  in  (10):  Condition  and  FjFjG15.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  96:  PDC  as  in  (10):  Condition  ^3  and  EEG15. X-aa;es  in  all  plots  are 
freqnencies  0  —  20Hz  (DC  and  freqnencies  aronnd  IHz  are  usnally  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  97;  PDC  as  in  (10):  Condition  and  FjEGlG.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  98:  PDC  as  in  (10):  Condition  ^3  and  FjEGlG.X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 


57 


1 


C4-A1 


C4-A1 


C4-A1 


C4-A1 


CJ 

LL 

^  0.5 

Q. 

Oz-A 

o 

cn 

itfliHB 

0 

^ - 

0 

— 

'  0 

10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


C3-A2  C3-A2  C3-A2  C3-A2 


^  or  :  :  :  '  0^^  ■  ■  '  0^^:  ^ 

10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 

Fz-AI  Fz-A1  Fz-AI  Fz-AI 


1| - 1  1 


10  20  30  40  °  10  20  30  40 

Pz-A1  Pz-AI 


Pz-AI 


CJ 

< 

0.5 

lllll  1  IIIIIIIIIIIMtt^^ 

0 

'  . rTTn|  ^  I . I 

10  20  30  40  10  20  30  40 

Pz-AI 


1| - 1  1| - 1  1| - 1  1 


10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 

Oz-A2  Oz-A2  Oz-A2  Oz-A2  EEG17;1;PDC 


C\J 

< 

cJ.  0  5 

^  0.5 

N  0-5 

o 

0 

LL 

0 

_ 

Q. 

0 

10  20  30  40  10  20  30  40  10  20  30  40  10  20  30  40 


Figure  99;  PDC  as  in  (10):  Condition  and  Yj'EGVJ .X-axes  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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Figure  100:  PDC  as  in  (10):  Condition  ^3  and  EEG17.X-aa;es  in  all  plots  are 
frequencies  0  —  20Hz  (DC  and  frequencies  around  IHz  are  usually  filtered  by  EMD 
thus  no  activity  there  analyzed). 
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4  Future  work 


4.1  Statistical  analysis  of  synchrony  patterns 

Confidence  levels  of  synchrony  “patterns”  in  around  20Hz  and  around  70Hz  to  be 

analyzed. 

4.2  Automatic  classification/recognition  based  on  synchrony 
patterns 

4.2.1  Two  classes  of  no-nap  and  40  minute  nap  conditions 

4.2.2  Trial  of  clustering/classification  of  test  batteries  8  to  17 
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